function output = fn_mg_cost(delta,R1,D0,chi,s,phi,dist_d,d_min,d_max,kappa)

kappaprime = kappa(1)*exp(kappa(2)*fn_tax(delta,R1,D0,chi,s,phi,dist_d,d_min,d_max));

output = -(1 + kappaprime);

end